Infrared conductivity of a one-dimensional charge- ordered state: 

quantum lattice effects 

C.A. Perroni, V. Cataudella, G. De Filippis, G. Iadonisi, V. Marigliano Ramaglia, and F. 

Ventriglia 

Coherentia-INFM and Dipartimento di Scienze Fisiche, 
Universitd degli Studi di Napoli "Federico II", 
Complesso Universitario Monte Sant'Angelo, 
Via Cintia, 1-80126 Napoli, Italy 

(February 2, 2008) 
The optical properties of the charge-ordering (CO) phase of the one- 
dimensional (1-D) half-filled spinless Holstein model are derived at zero tem- 
perature within a well-known variational approach improved including second- 
order lattice fluctuations. Within the CO phase, the static lattice distortions 
give rise to the optical interband gap, that broadens as the strength of the 
electron-phonon (el — ph) interaction increases. The lattice fluctuation ef- 
fects induce a long subgap tail in the infrared conductivity and a wide band 
above the gap energy. The first term is due to the multi-phonon emission 
by the charge carriers, the second to the interband transitions accompanied 
by the multi-phonon scattering. The results show a good agreement with 
experimental spectra. 
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In the last years there has been a renewed interest in charge density wave (CDW) 
materials. 1 The transition to a CO phase is common to a wide range of compounds, 2 in- 
cluding many quasi-l-D materials such as organic conjugated polymers, charge transfer salts, 
molybdenum bronzes and MX chains. 1-5 These materials undergo a Peierls instability driven 
by the el —ph interaction in the half-filled band case. For most quasi-l-D systems, the lattice 
zero-point motion is comparable to the Peierls lattice distortion, so quantum lattice fluctu- 
ations must be taken into account to satisfactorily describe spectral, transport and optical 
properties. 6-8 In particular the study of the optical absorption of CO materials represents 
a very useful tool to extract the gap energy and, in general, to investigate the properties of 
the ordered state. 1,9-11 In this framework, the experimental measurements have pointed out 
that the lattice fluctuation effects can remove the inverse square-root singularity expected 
for the case of a static distorted lattice 12 affecting profoundly the conductivity spectra. 6 ' 13 ' 14 
Actually these effects give rise to the subgap optical absorption seen in these materials where 
a significant tail below the maximum of the interband transition term is measured. Moreover 
the optical absorption above the interband optical gap band also presents deviations from 
the behavior obtained within the mean-field approach of the static lattice. 13 

The challenge of understanding the effect of quantum lattice fluctuations on the Peierls 
dimerization and the absorption spectra has determined an intense study of the Holstein 
model 15 that is a typical el — ph coupling model developing a CO state at half-filling. 
Actually the Holstein Hamiltonian has been investigated by using various techniques, such as 
Monte-Carlo simulations, 16 ' 17 renormalization-group analysis, 18 ' 19 variational method, 20-22 
density- matrix renormalization group 23 and exact diagonalization. 24 These studies reveal 
that, in the spinless case, there is a quantum phase transition from a Luttinger liquid 
(metallic) phase to an insulating phase with CDW long-range order. Because of limited 
system sizes in numerical approaches, except for the antiadiabatic regime, the behavior 
of the conductivity spectra is not well determined, so the extraction of the gap from the 
optical data is not precise. 24 However, through exact diagonalization methods, the spectral 
weight of the conductivity can be deduced showing the onset of the infrared absorption for 
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lattices of a few sites. Recently, an analytical variational approach, 21 ' 22 valid in the weak 
el — ph coupling regime, has been developed to study the phase diagram and the optical 
conductivity of the Holstein model. The employed approximation has different effects on the 
phase diagram and the optical properties. In fact, the maximum in the optical spectra is not 
directly corresponding to the gap calculated within the variational approach. Actually the 
peak position of the optical conductivity is higher than the gap with an energy separation 
of the order of the electron transfer integral in the adiabatic limit. Therefore, experimental 
findings of a tail in the optical spectra below the gap energy (corresponding to the maximum 
of experimental conductivities) does not find a clear explanation within previous approaches. 

In this paper we employ the variational scheme proposed for the ID half-filled spinless 
Holstein model by H. Zheng, D. Feinberg and M. Avignon 20 (ZFA) in order to calculate 
the spectral properties and the infrared response. We note that the ZFA method improves 
the mean- field solution showing a partial agreement with other numerical works 16 ' 17 ' 23 ' 24 . 
Actually this approach is able to introduce lattice fluctuations on the mean-field Peierls 
solution since it takes into account the nonadiabatic polaron formation. Indeed, due to the 
polaronic effect, the lattice deformation is allowed to follow instantaneously the electrons. 
The calculated conductivity spectra are characterized by a transfer of spectral weight from 
low to high energies and by a broadening of the optical interband gap, with increasing the 
el — ph coupling. The effect of the quantum lattice fluctuations is able to determine in 
the infrared conductivity a subgap absorption term near the phonon energy and a wide 
band above the gap energy. The first contribution is due to the multi-phonon emission 
by the charge carriers, the second to the interband transitions accompanied by the multi- 
phonon scattering. The inclusion of lattice fluctuation effects beyond the ZFA approach is 
able to smooth the inverse square-root singularity of the ZFA and mean-field conductivity. 
Moreover these effects strengthen the features of the conductivity below and above the 
gap already found in the ZFA approach inducing an actual subgap tail. Therefore lattice 
fluctuations influence the optical absorption at the low-energy subgap, at the gap and at 
the high-energy scale above the interband optical gap term. These features are found in the 
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measured spectra. 13 ' 14 

In section I the model and the ZFA variational approach are briefly reviewed; in section 
II the spectral properties are deduced in order to characterize the quasi-particle gap of 
the CO phase; in section III the infrared spectra obtained within the ZFA approach are 
discussed; finally in section IV the effects of lattice fluctuations beyond the ZFA approach 
are analyzed. 

I. MODEL AND ZFA VARIATIONAL APPROACH 

In this paper we deal with the ID spinless Holstein model 15 at half-filling. 
The Holstein Hamiltonian is 

H = -t C i C j + ^0 J2 a l a i + 9^0 J2 C \°i ( a i + a i) - A* S C i Ci ' W 

<i,j> i i i 

Here t is the electron transfer integral between nearest neighbor (nn) sites < i,j >, cj (q) 
creates (destroys) an electron at the i-th site and \x is the chemical potential. In the second 
term of eq.(l) a\ (a*) is the creation (annihilation) phonon operator at the site i, uq denotes 
the frequency of the optical phonon mode and the parameter g represents the coupling 
constant between electrons and phonons. The dimensionless parameter A 

A = (2) 

indicating the ratio between the small polaron binding energy and the energy gain of an 
itinerant electron on a rigid lattice, is useful to measure the strength of the el—ph interaction 
in the adiabatic regime. We consider spinless electrons since they, even if at a very rough 
level, mimic the action of an on-site Coulomb repulsion preventing the formation of local 
pairs. Actually, for one dimensional systems in the limit of infinite local repulsion U, the 
charge sector of the Hubbard model maps onto a spinless model, therefore the spinless 
Holstein model can be considered as a reliable model for typical one-dimensional systems. 
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The hopping of electrons is supposed to take place between the equivalent nn sites of a 
ID lattice separated by the distance a. The units are such that the Planck constant k — 1 
and the lattice parameter a—1. 

As stressed in the original ZFA paper, 20 the starting point of the approach is the con- 
sideration that the strong coupling and infinite phonon frequency limit of the model are 
described by polarons. The ZFA approach extends the polaron formation to intermediate 
regimes recovering the mean-field solution of the zero phonon frequency limit. Following the 
ZFA variational scheme, three successive canonical transformations are performed in order 
to treat the electron-phonon interaction variationally and to introduce the charge-ordering 
solution. 

The first transformation is the variational Lang-Firsov unitary one 25 



Ui = exp 



(3) 



where / and Aj are variational parameters. The quantity / controls the degree of the an- 
tiadiabatic polaronic effect since the lattice deformation is allowed to follow instantaneously 
the electrons. Moreover Aj denotes a displacement field describing lattice distortions due 
to the average electron motion. At half-filling the charge-ordered solution is obtained by 
assuming 

A< = A + A co e iQR % (4) 

where A represents the lattice distortion unaffected by the instantaneous position of electrons 
and Aco the additional local lattice distortion due to the Peierls dimerization with Q = n. 
The second transformation is 



U 2 = exp 



a ^2 ( a ] a i — a j a j 



(5) 
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where the variational parameter a determines a phonon frequency renormalization. 
The transformed Hamiltonian H = U^U^ 1 HU\U2 is 
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a = -t e x \ x A c 3 + E + Luj o sinh2 ( 2 «) + 9 2 u E A * 2 

<i,j> i i 

+u>o sinh (2a) cosh (2a) E + ciiCLij — g^oe 2a E A « + a t) 

i i 

+gu (1 - /) e 2a E cjci (a,i + af) + E cjci (^ - //) , (6) 

i i 

where we have introduced the phonon operator 

X, = exp [gfe- 2a (a t - a})} , (7) 

the renormalized phonon frequency u> = cosh(4a) , the number of lattice sites L and the 
quantity rji 

Vi = <?W (/ - 2) + 2g 2 cu (/ - 1) Aj. (8) 

In the ZFA approach the energy is deduced introducing a test Hamiltonian characterized 
by non interacting electron and phonon degrees of freedom such that (H — H test ) t = 0, 
where <> t indicates the mean value obtained by using the ground state of H test . The test 
Hamiltonian is given by 

Htest = -t e ff E c l c i + E a l 0i + Loj v sinn2 ( 2a ) + L 9 2 ^o (a 2 + A^ ) 

<i,j> i 

-2g 2 u A co (1 - /) E e iQRi c\ Ci - ^ £ c ! c » ( 9 ) 

i i 

where the subsidiary chemical potential /xo is 

/x = /i - <?W (/ - 2) + 2^^ (/ - 1) A. (10) 
The quantity t e ff = te~ s denotes the effective transfer integral, where the quantity 

S = g 2 fe-^ (11) 

controls the band renormalization due to the nonadiabatic polaron effect. The electronic 
part of the test Hamiltonian is diagonalized by a third canonical Bogoliubov transformation 20 
yielding 



6 



Htest = - fMi) d\d k + (d } - Mo) PkPk + 51 4 a </ + 

fceTVZ keNZ q 

+Lu sinh 2 (2a) + Lg 2 u (A 2 + A^ D ) , (12) 



where d\ (dk) creates (destroys) a quasi-particle in the upper band £^ = ^Je 2 + E 2 , p|. (p^) 
creates (destroys) a quasi-particle in the lower band = — ^Je 2 k + E 2 and e fe is the polaronic 
band. We note that NZ indicates the New (Brillouin) Zone defined by the condition e k < 0. 
In the CO phase a gap opens between the upper and lower bands in the quasi-particle 
spectrum and it is twice the quantity 

E = 2g 2 u A co (1 - /) • (13) 

Within the variational approach the kinetic energy mean value E kin is 

E km =<f>=- f°_ deg{e)-^==, (14) 
J-w Ve 2 + E 2 

where W = 2t e ff is the renormalized band half-width, g(e) the ID density of states, and 
the electron order parameter m e is given by 



-) T,e&*(c\ci) = E f° de ^ . (15) 
LJ Y J-w Ve 2 + E 2 V ; 



The CO phase is characterized by the order parameter different from zero. 

In Fig.l we report the phase diagram 20 derived within the ZFA approach in the ther- 
modynamic limit. The ordered state is separated with a transition line by the A phase that 
represents the disordered phase (A C o = 0). Previous studies 17,21,23 ' 24 have pointed out that 
the normal state has the properties of a Luttinger liquid. This has been verified also by using 
the ZFA wave-function and making a finite-size scaling analysis. 24 In the inset of Fig.l, there 
is the comparison between the transition lines calculated in a mean-field approach (white 
squares), ZFA approach (white circles) and DMRG approach 23 (black squares). Here mean- 
field approach means that we are neglecting the effect of polaron formation (/ = 0, a — 0). 
In the ranges of parameters relevant for quasi one-dimensional materials, the ZFA approach 
improves the mean-field solution since the CO phase is stable for larger values of the el—ph 



couplings. However the lattice fluctuation effects introduced by this approach are not suffi- 
cient to obtain transition lines comparable with those of DMRG approach. This indicates 
that quantum fluctuation effects beyond the ZFA scheme are not negligible and can play a 
role also in the calculation of the dynamical properties of the model. 



II. SPECTRAL PROPERTIES WITHIN ZFA APPROACH 



In this section we calculate the spectral properties within the ZFA approach. They are 
discussed in order to characterize the gap in the quasi-particle spectrum. 

The electron retarded Green's function can be disentangled into electronic and phononic 
terms 26 by using the test Hamiltonian (9), hence 



G ret {k,t) = e~ s G { ^{k,t) + e- s [e^{Se^ t )-l]- £ G^(h,t) + 



kxeNZ 



eq> (Se***) - l] - ]T G ( -J(h,t). 



(16) 



k^NZ 



In equation (16) G^(k,t) is 



Gret{k>t) — u k^ret {k,t) + V k G ( r J (k,t) , 



,2 n (+) 



2^(-), 



(17) 



where G [ ^(k : t) and G { ~}{k,t) are the Green's functions associated to the quasi-particles of 
the upper and the lower bands, respectively, with u\ given by 

1 



1 + 



€k 



?2 



+ E 2 



(18) 



and v\ by 



v k = - 2 



?2 



+ E 2 



(19) 



One obtains the spectral function 



A(k, u) = -2ZG ret (k, u) = 2ne~ s [u 2 k 5 (u - £ ( k +) ) + v 2 k 5 (u - £ ( k +) ) 



+ 



S 



2ne s ^2 — r [H(w — ncu ) + H(—cu — ucuq)} , 



n=0 



nl 



(20) 



where the function H{uj) is 



H{u) 



g (Vuji - E 2 ) 



6(u-E)6 (VE 2 + W 2 - u) , 



(21) 




with 6{x) Heaviside function. Two physically distinct terms 26 appear in eq.(20): the coherent 
and incoherent one. The first derives from the coherent motion of charge carriers and their 
surrounding phonon cloud. In the normal phase it represents the purely polaronic band 
contribution and shows a delta behavior. In the CO phase, this term is equal to the result 
of the mean-field approach when one neglects the renormalization of the upper and lower 
band due to the polaron effect. The incoherent term in eq.(20) is due to processes changing 
the number of phonons during the hopping of the charges, and provides a contribution 
spreading over a wide energy range. 

In Fig. 2 we report the renormalized density of states N{u) calculated in the ZFA ap- 
proach (solid line) and mean-field (dashed line) at a fixed value of the el — ph coupling and 
t = 5uq. In the CO phase, a gap opens in the quasi-particle spectrum and it is larger for the 
mean-field solution than the ZFA one. We note that, at the energies corresponding to the 
gap, the inverse square-root singularity occurs for both approaches. The other sharp feature 
in the density of states derived in the ZFA approach is due to one-phonon processes in the 
upper and lower bands that are relevant in the intermediate el — ph coupling regime. 



In this section we focus our attention on the optical properties within the ZFA approach. 
Since we are primarily interested to the absorption spectra in the CO phase, we evaluate 
the conductivity for the frequency uj different from zero. 

In a regime of linear response the real part of the conductivity is given by the current- 
current correlation function 



III. OPTICAL PROPERTIES WITHIN ZFA APPROACH 




(22) 
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where (5 is the inverse of the temperature and j the current operator. Performing the two 
canonical transformations (3,5) and making the decoupling 26 of the correlation function in 
the electron and phonon terms through the introduction of H test (9), we get 



<J f Wj(0)> = E E( 5 • W, i\ S, 6', t)A(i, i', 5, 6', t), 

i,S i',8' 

where the function A (i, i', 5, 5', t) denotes the electron correlation function 



A (i, i', 5, 5', t) = {<i(t)c i+ s(t)4> +S 'Ci')t 



and the function $ (i, i', 5, 5', t) the phonon correlation function 



(23) 



(24) 



*(i,i',6,6',t) = (X}(t)X l+s (t)X} /+s ,X t ,) t . 



In order to simplify the analysis of our results, we separate $ into two contributions 



(25) 



$ (i, i', 5, 5>, t) = [(XjX t+5& ) t \ |$ (i, i', 5, 5', t) - \{x}x i+ sa) t 

$ (i,i',S, 5',t)-e- 2S 



= e~ 2S + 



(26) 



where S is given by eq.(ll). 

Considering eq.(26), the conductivity can be expressed as a sum of two terms 26 

Ua{u) = ^a {coh \u) + 9&7 (inooh) (w). 



(27) 



As in the spectral properties, the appearance of two physically distinct contributions, the 
coherent and incoherent one, occurs. Actually the first term ^R,a^ coh ^ is due to the charge 
transfer affected by the interactions with the lattice but not accompanied by processes 
changing the number of phonons. On the other hand, the incoherent term $la( mcoh ^ derives 
from inelastic scattering processes of emission and absorption of phonons. Both terms of the 
conductivity can be expressed in terms of the Green's function G^(k,t) given in eq.(17). 
The coherent conductivity is derived as 

fA7ie 2 t r 



to 



2-2S 



LO 
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with f (+) = Ve 2 + E 2 . We note that this term gives contribution to the conductivity only 
in the CO phase since it directly depends on the semi-gap E. Furthermore it is equal to 
the mean-field conductivity 12 when the renormalization of the upper and lower bands due 
to the polaron effect is neglected. 

The incoherent term of the conductivity can be divided into two components: 

The quantity $lai ncoh \uj) is due to the multi-phonon emission by the charge carriers in the 
lower band £^ that does not change the electron momentum. This first term reads 



(incoh) / \ 27Te t _ 2 g 
UJ 



\ S cos ( k )( v l -«fc) 



1 2 



E^^-^o), (30) 
n! 



keNZ J n=l 

that, making the envelope of the delta functions (procedure exact in the limit u — > 0), 
becomes 

^'M^^rlrS)^-^' (31) 

where E kin is the mean value of the kinetic energy equal to eq.(14) and T(x) is the gamma 
function. 

In eq.(29), Kcr^"^^) takes into account the interband transitions accompanied by 
multi-phonon scattering. This second term is given by 

W L k u k 2 eNZ n=l H - 

(32) 

that, enveloping the delta functions, can be transformed as 

» (inooft), x /27re 2 t 2 \ _ 2S f o f o (2gp f E 2 \ 

>) = (^-J e U dei U de2 wTV) 9{ei)9{e2) { 1+ ww 

x9(u - u) - 2E) (33) 



where 7\ = min - cj - E) 2 - E 2 ) , T 2 = min (w, \j {u — u Q — ^) 2 — E 2 ^ , ^ +) 



is defined as 
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with i = 1, 2, and y = (w - - Q )M)- 
We have checked that the sum rule 

/ dutta(u) = --e 2 E km (35) 

is verified by the calculated conductivity spectra in the CO phase, where E kin is given by 
eq.(14). 

In Fig. 3a we report the different contributions to the conductivity spectrum derived 
within the ZFA approach for t = 5uj and g = 2.4: the coherent term (dashed line), the 
incoherent term due to multi-phonon emissions (dotted line) and the incoherent term in 
correspondence with interband transitions (dash-dotted line). For this value of the el — 
ph coupling, the gap is larger than the phonon frequency, that represents the absorption 
threshold of the first incoherent term. The incoherent term at Uo reaches the largest values in 
the intermediate coupling regime where the nonadiabatic polaron formation is more effective. 
This contribution is still present for large el—ph couplings giving rise to a subgap absorption 
band. The second term of the incoherent conductivity spreads for a wide range of frequencies 
above the energy gap following the coherent interband absorption band. 

In Fig. 3b the conductivity spectra within ZFA (solid line) and mean-field (dashed line) 
approach at t — 5u and g = 2.4 are compared. Not only the optical gap within ZFA 
is smaller than mean-field one, but the two incoherent terms due to lattice fluctuation 
effects are able to provide a non negligible contribution below and above the gap (the arrow 
in figure indicates this last contribution). This suggests that a better treatment of the 
lattice fluctuations allows to capture features of the optical conductivities that are found in 
experimental spectra. 13 In the next section we will see that lattice fluctuation effects beyond 
the ZFA approach give rise in the conductivity to an actual subgap tail in good agreement 
with experimental data. 

With rising the el — ph coupling, a transfer of spectral weight from low to high energies 
takes place and the optical gap broadens. The optical response within ZFA approach is 
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strongly dependent on adiabaticity ratio since the incoherent terms acquire increasing spec- 
tral weight compared with that of the coherent interband term as the ratio t/cu decreases. 
Actually the spectral weight can be measured through the u;-integrated function 



In this section we deal with quantum lattice fluctuation effects beyond the ZFA ap- 
proach. We first determine the scattering rate of the quasi-particles of the upper and lower 
bands. Next we will analyze the effect of the self-energy insertions on the infrared conductiv- 
ity that is the main aim of this paper. This is performed following the lines of our previous 
works. 26 ' 27 

Now we consider the actual transformed Hamiltonian of the system in eq.(6) as a per- 
turbation of the test Hamiltonian in eq.(9). At the second order of the perturbation theory 
the retarded self-energy (k, u) can be derived. 28-32 If only the dominant autocorrela- 
tion terms are retained, the self-energy is local and provides the scattering rates of the 
quasi-particles of the upper and lower bands 




(36) 



whose value S m in the limit of infinite frequency is given by eq.(35). 



IV. FLUCTUATIONS BEYOND ZFA APPROACH 




(37) 



These two quantities turn out to be equal, so we have only one scattering rate 



r(fc) = r + (fc) = r_(fc). 



(38) 



The scattering rate can be decomposed as 



r(fc) — r(£ fc ) — Ti-phoniik ) + ^ multi-phon{ik )) 



(39) 



where Ti_ phon is the contribution due to single phonon processes only 
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1.2-2SJ2 f2-4a 



(+)^ 



-.2, ,2„4ai 



(40) 



Tmuiti-phon represents the scattering rate by multiphonon processes 

+°° (Or, 2 f2„-4a\ l 

1 multi-phon\Kk ) — Zt e 2^ 7j 01ASfc 

1=2 L 

In the previous equations the function gij(£,^) reads 



(41) 



9i,i(d +) ) = n F (£F + too) ^(^ ; + too) 



*(+) 



+ 



l-^(d +) -to ) ^r-too 



(42) 



and 5f 2 (£ 



(+)> 



92(& +h 



n F (d +) +^o) 5(^+^0)+ l-n F (# 



r(+) 



(+) 



(43) 



where B(x) = 2n [H{x) +H(— x)], with H(x) given in eq.(21). In the CO phase the scattering 
rate has a gap due to the dimerization and the process of phonon spontaneous emission by 
the quasi-particles. 27 

The role of the scattering rate is important not only to improve the approximations of 
calculation of the spectral properties, but also the optical properties. Through the scattering 
rate, we can consider the new Green's function G^(k,t) 



G&\k,t) = ulGl + J(k,t)+vlGU(k,t), 



,2A(+) 



(44) 



that substitutes the function G^{k,t) given in eq.(17). In eq.(44) the Green' functions 



gQ(M) are 



G%(k, t) = -iB(t) exp (-itPt) exp {-tT(k)/2) , 



(45) 



with v standing for + or — . Therefore it is possible to derive a new spectral function 
and density of states that include lattice fluctuation effects beyond the ZFA approach. A 
pseudogap as precursor of the actual gap at stronger el —ph couplings is found in the density 
of states. 27 
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The inclusion of the scattering rate is able to affect the features of the infrared absorption. 
As in the ZFA approach, the conductivity is decomposed into a coherent and an incoherent 
term: 



Ua Fluc M = Ka^UlM + H*%£\u>). (46) 



The coherent conductivity is given by 



= (—) ^ E f- de[n F {^ - U ) -n F (^))]C(-'-)(6,a;)M6)A(-'-)(6), 



(47) 

where C^ V2 \e,uj) is 

C^\e,u) = , F(e) - 2, (48) 

v ; r 2 (e) + (^ -e vi) + ^f 

h{e) = g{e) (l — -Jz—), with gr(e) bare density of states, and the function A^ Ul ' V2 \e) is 
expressed by 

A(+'+)(c)=A(-'-)(c) = ^ 5 (49) 



with 



A (+ '- ) (c)=A(-'+)(c) = F ^. (50) 



The latter term of the conductivity becomes 



= ( — ) e" 2 * 1 E f - de f - de l9 (e)9(e 1 )R^\e,e 1 ,uj), (51) 

\ U J V1,V2 J ~ W J ~ W 



where the function R^ Ul ' V2 \e, ei,u) is given by 

+°° (2g 2 f 2 e- 4a ^ 1 
i=i 

C^' U2 \e,e u x) is 



i?(-'-)(e, ei ,c) = E 1 ^ ^ [^'"Wi^) + H}™\e,e 1 ,u>)] , (52) 



C («) (e £] x) = 1 F(6)+r( Cl )] 

2 [r ( e ) + r( ei )]V4+(^)-^ 2) + a:) 2 



(53) 
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and 



(54) 



In eq.(51) the functions J; (e, ei,w) 

^(Mlw) = C^\e,e 1 ,u + lu )[n F {^ - loo -co)- n F (^ - lu )]n F (& 2) ) (55) 
H^ 2 \e, Cl , w ) = C^fo Cl , w - to )[^(^ 2) + tob - w) - + fob)] x 



describe phonon scattering processes. 

In Fig. 4 we show the conductivity at t — 10u and A = 0.45: the solid line represents 
the quantity obtained adding fluctuations over ZFA, the dashed line that derived within 
ZFA and dotted line that of the mean-field approach. The values of the parameters have 
been chosen such that they are appropriate for quasi ID inorganic metals. First we observe 
that the inverse square-root singularity obtained in the mean-field and ZFA approach is 
reduced. Then, due to quantum lattice fluctuations, the new conductivity not only shows a 
long subgap tail but presents anomalies also over the gap when compared with the mean-field 
spectrum. Therefore it strengthens the tendencies already shown by the ZFA conductivity 
and it is in good agreement with the experimental spectra. 6 ' 13 ' 14 

As stressed in the introduction, an analytical variational approach, 21,22 valid in the weak 
el—ph coupling regime and based on a similar procedure of calculation, has been developed to 
study the optical conductivity. In order to emphasize the different physical results between 
the present work and the previous one, in Fig.5 we compare the conductivity ^R,<Jpi uct (solid 
line) with the corresponding quantity (dashed line) obtained within the preceding approach 22 
for t = 10uo and A = 1. This last conductivity is extrapolated from the weak coupling 
limit and, as already noted by the authors of the work, it shows a maximum that does not 
coincide with their energy gap (see dashed arrow in the figure) . This feature is not consistent 



(56) 
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with experimental spectra 13 that, at T = 0, show the maximum of the infrared optical 
conductivity generally near to the gap determined by other experimental measurements 
such as the low-temperature resistivity, the temperature dependence of the susceptibility and 
neutron and Raman scattering. 14 As clearly shown in Fig.5, our approach is not affected by 
this problem (the gap in our approach is indicated by the solid arrow). While in the previous 
work, 22 the presence of a tail below the peak should be due to usual interband transitions, 
in our approach the peak in the conductivity corresponds to the gap in the quasi-particle 
spectrum and is well above a real subgap tail determined by lattice fluctuations. 

In the inset of Fig.5 we report the ratio (solid line) between S(u), the spectral weight 
of the conductivity calculated in this section, and S m , the same quantity in the limit of 
infinite frequency. This ratio is compared with that (dashed line) obtained by means of 
exact diagonalizations of the Hamiltonian for a system of six sites. 24 Our approach captures 
the correct onset of the optical response, even if it presents sharper features. This could 
be due to the fact that we are not considering a small cluster but we are performing the 
thermodynamic limit, and, furthermore, that higher order self-energy insertions should be 
included in order to improve the approach. 

The optical response shows similar features at lower values of the adiabaticity ratio 
t/ujo- In the regime of near electronic and phononic energy scales the effects due to lattice 
fluctuations are enhanced and the contributions below and above the gap subtract a larger 
spectral weight to the interband gap term. This regime is typically important for inorganic 
linear chain compounds, for example, the compound TTF — TCNQ, that is a narrow-band 
one-dimensional metal with a relatively strong el — ph coupling. 33 

V. CONCLUSIONS 

We have discussed the optical properties of the half-filled spinless Holstein model in 
the ID case at zero temperature within the ZFA approach. We have observed that, with 
increasing the el — ph coupling, the ordered phase affects the conductivity spectra inducing 
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a transfer of spectral weight from low to high energies and a broadening of the optical 
interband gap. The quantum lattice fluctuations considered in the ZFA approach are able 
to affect the optical response of the system that shows bands below and above the gap. 
When fluctuation effects beyond the ZFA approach are included, the optical conductivities 
are profoundly changed since they are characterized by a subgap tail and a wide band above 
the interband optical gap term that is smoothed when compared with the mean-field result. 
Therefore the inelastic scattering processes influence the low-frequency and high-frequency 
features of the conductivity in agreement with experimental spectra. 13 ' 14 Our results make 
clear that a treatment of the lattice fluctuations beyond the ZFA approach is required to 
obtain a consistent agreement with experimental data. 

In this paper lattice fluctuation effects beyond the ZFA approach are included calcu- 
lating perturbatively the scattering rate of the polarons that form the ordered state. A 
second-order perturbation calculation on the ZFA solution reproduces the integrated spec- 
tral weights obtained by exact numerical approaches suggesting that the employed approach 
can capture the infrared response of the model. Finally we note that our approach is valid 
in the infrared range of frequencies where the interband absorption typically takes place. 
Thus it is not able to reveal the structures attributed to collective excitation modes arising 
from the CDW condensate. 1 ' 9,10 



FIGURE CAPTIONS 

Fl The phase diagram of the system. The line separates the CO from the A phase, that 
represents the disordered normal state. In the inset, the transition line between the 
CO and the disordered A phase calculated in the mean-field (white squares), ZFA 
(white circles) and DMRG approach (black squares). 

F2 The renormalized density of states derived in the ZFA (solid line) and mean-field 
approach (dashed line) at t — 5u and g = 2.2 as a function of the energy (in units of 

Uo). 
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F3 (a) The conductivity at t — 5uo and g = 2.4 as a function of the frequency (in units 
of u ) decomposed into its different components: coherent term (dashed line), first 
incoherent term (dotted line) and second incoherent term (dash-dotted line). 

(b) The conductivity spectra derived in the ZFA (solid line) and mean-field approach 
(dashed line) at t — 5u and g = 2.4 as a function of the frequency (in units of uq). 

The conductivities are expressed in units of with m = ^. 

F4 The conductivity (in units of with m = ^) obtained including fluctuations beyond 
the ZFA (solid line), derived in the ZFA (dashed line) and in the mean-field approach 
(dotted line) up to 14 uo at t = IOuq and A = 0.45. 

F5 The conductivity of the present approach (solid line) in comparison with that (dashed 
line) calculated in a previous paper 22 at t — 10u and A = 1.0. The solid arrow 
indicates the gap in our approach and the dashed arrow that obtained in a previous 
work. 22 In the inset the ratio between the spectral weight S(u) and S m calculated 
in the present approach (solid line) compared with the same quantity (dashed line) 
derived from exact numerical diagonalizations. 24 
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